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Abstract 

> 

Using the equation of state of asymmetric nuclear matter that has been re- 
cently constrained by the isospin diffusion data from intermediate-energy heavy 
ion collisions, we have studied the transition density and pressure at the inner 
edge of neutron star crusts, and they are found to be 0.040 fm~ 3 < pt < 0.065 
fm" 3 and 0.01 MeV/fm 3 < P t < 0.26 MeV/fm 3 , respectively, in both the dy- 
namical and thermodynamical approaches. We have further found that the 
widely used parabolic approximation to the equation of state of asymmetric 
nuclear matter gives significantly higher values of core-crust transition density 
^ . and pressure, especially for stiff symmetry energies. With these newly deter- 

mined transition density and pressure, we have obtained an improved relation 
between the mass and radius of neutron stars based on the observed minimum 
crustal fraction of the total moment of inertia for Vela pulsar. 
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1 Introduction 

Exploring the properties of neutron stars, which are among the most mysterious 
objects in the universe, allows us to test our knowledge of matter under extreme 
conditions. Theoretical studies have shown that neutron stars are expected to have a 
liquid core surrounded by a solid inner crust [1] , which extends outward to the neutron 
drip-out region. While the neutron drip-out density p out is relatively well determined 
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to be about 4 x 10 11 g/cm 3 or 0.00024 fm -3 [2], the transition density pt at the inner 
edge is still largely uncertain mainly because of our very limited knowledge on the 
nuclear equation of state (EOS), especially the density dependence of the symmetry 
energy (E sym (p)) of neutron-rich nuclear matter (3JH]. These uncertainties have ham- 
pered our understanding of many important properties of neutron stars [31 HI [5] and 
related astrophysical observations [31 [H El El El El HOI HH [12l [131 H3] . Recently, signif- 
icant progress has been made in constraining the EOS of neutron-rich nuclear matter 
using terrestrial laboratory experiments (See Ref. [15] for the most recent review). In 
particular, the analysis of the isospin-diffusion data [T6l [T71 [T8] in heavy-ion collisions 
has constrained tightly the E sym (p) in exactly the same sub-saturation density region 
around the expected inner edge of neutron star crust. The extracted slope parameter 
L = 3p ag g m |p=p m the density dependence of the nuclear symmetry energy was 



found to be 86 ± 25 MeV [19J, which has further been confirmed by a more recent 

analysis [20J. With this constrained nuclear symmetry energy and the corresponding 
EOS of asymmetric nuclear matter, we have obtained an improved determination of 
the values for the transition density and pressure at the inner edge of neutron star 
crusts. This has led us to obtain significantly different values for the radius of the Vela 
pulsar from those estimated previously. Also, we have found that the widely used 
parabolic approximation (PA) to the EOS of asymmetric nuclear matter gives much 
larger values for the transition density and pressure, especially for stiff symmetry 
energies. 

2 Dynamical and Thermodynamical approaches to 
the stability of npe matter 

The inner edge of neutron star crusts corresponds to the phase transition from the 
homogeneous matter at high densities to the inhomogeneous matter at low densities. 
In principle, the inner edge can be located by comparing in detail relevant properties of 
the nonuniform solid crust and the uniform liquid core mainly consisting of neutrons, 
protons and electrons (npe matter). However, this is practically very difficult since the 
inner crust may contain nuclei having very complicated geometries, usually known as 
the 'nuclear pasta' [3, [12j [211 1221 [23] . Furthermore, the core-crust transition is thought 
to be a very weak first-order phase transition and model calculations lead to a very 
small density discontinuities at the transition [21 EH ESI [25] ■ In practice, therefore, 
a good approximation is to search for the density at which the uniform liquid first 
becomes unstable against small amplitude density fluctuations with clusterization. 
This approximation has been shown to produce very small error for the actual core- 
crust transition density and it would yield the exact transition density for a second- 
order phase transition 0, EH [251 EE]. Several such methods including the dynamical 
method El El EH ETJ [28] , the thermodynamical method (HEH1I30] and the random 
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phase approximation (RPA) [3U [26] have been applied extensively in the literature. 
Our study here was based on the dynamical and thermodynamical approaches. 

In the dynamical approach, the stability condition for a homogeneous npe matter 
against small periodic density perturbations can be well approximated by P, [7J, El El 

47re 2 

VVO) = v + (3k 2 + > 0, (1) 

where k is the wavevector of the spatially periodic density perturbations and 
v = (dp n /dp p ) 2 k2 = Aire 2 

° dp p dp n /dp n ' TF dp e /dpe 



(3 = D pp + 2D np C + D nn C 2 , C = -f^, (3) 

dUn/dpn 



with pi being the chemical potential of particle type i. The three terms in Eq. ([I]) 
represent, respectively, the contributions from the bulk nuclear matter, the density- 
gradient (surface) terms, and the Coulomb interaction. For the coefficients of density- 
gradient terms, we use the empirical values of D pp = D nn = D np = 132 MeV-fm 5 , 
which are consistent with those from the Skyrme-Hartree-Fock calculations [271 [32] . 
At k min = l(^f-) 1/2 - kj, F ] 1/2 , V dyn (k) has the minimal value of V dyn (k min ) = V + 
2(4vre 2 /3) 1 / 2 - fik 2 TF P El El [27J , and p t is then determined from V djn (k min ) = 0. 

In the thermodynamical approach, the stability conditions for the npe matter 
are H [2H E3] 

-© - - -(&>* 

In the above, P = Pj, + P e is the total pressure of the npe matter with P^ and P e 
being the contributions from baryons and electrons, respectively; v and q c are the 
volume and charge per baryon number; and p = p n — p p is the chemical potential of 
the npe matter. It can be shown that the first inequality in Eq. (j4j) is equivalent to 
requiring a positive bulk term Vo in Eq. (CQ) [32], while the second inequality in Eq. (jlj) 
is almost always satisfied. Therefore, the thermodynamical stability conditions are 
simply the limit of the dynamical one for k — > when the surface terms and the 
Coulomb interaction are neglected. 



3 Results for the transition density and pressure 

We have used in our study a momentum-dependent MDI interaction that is based 
on the modified finite-range Gogny effective interaction [31]. This interaction, which 
has been extensively studied in our previous work [15], is exactly the one used in 
analyzing the isospin diffusion data from heavy-ion reactions [T71 [T8] . 
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Figure 1: (Color online) The transition density pt as a function of the slope parameter L 
of the nuclear symmetry energy from the dynamical and thermodynamical approaches with 
and without the parabolic approximation in the MDI interaction. Taken from Ref. |32j . 

Shown in Fig. [1] is the transition density pt as a function of the slope parameter 
L of the symmetry energy from the MDI interaction. For comparisons, we have in- 
cluded results from both the dynamical and thermodynamical approaches with the full 
EOS and its parabolic approximation, i.e., E(p, 5) = E(p, 5 = 0) + E sym (p)S 2 + 0(5 4 ), 
where 5 = {p n ~Pp)/p is the isospin asymmetry, from the same MDI interaction. With 
the full MDI EOS, pt is seen to decrease almost linearly with increasing L in both 
approaches consistent with the results obtained in RPA pFlJ. Both dynamical and 
thermodynamical approaches give very similar results with the former having slightly 
smaller p t than the later (the difference is actually less than 0.01 fm -3 ), and this is due 
to the fact that the density-gradient and Coulomb terms in the dynamical approach 
make the system more stable and thus lower the transition density. The small dif- 
ference between the two approaches implies that the effects of density-gradient terms 
and Coulomb term are, however, unimportant in determining p t . On the other hand, 
significantly larger transition densities are obtained in the parabolic approximation, 
including the predictions by Kubis using the MDI EOS in the thermodynamical ap- 
proach [29], especially for stiffer symmetry energies (larger L values). The large error 
introduced by the parabolic approach is understandable since the /9-stable npe matter 
is usually highly neutron-rich and the contribution from the higher-order terms in 5 
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is thus appreciable. This is especially the case for the stiffer symmetry energy which 
generally leads to a more neutron-rich npe matter at subsaturation densities. Fur- 
thermore, because of the energy curvatures involved in the stability conditions, larger 
factors are multiplied to the contributions from higher-order terms in the EOS than 
that multiplied to the quadratic term. These features agree with the early finding 
[35] that the transition density p t is very sensitive to the fine details of the nuclear 
EOS. According to results from the more complete and realistic dynamical approach, 
the constrained L limits the transition density to 0.040 fm~ 3 < p t < 0.065 fm -3 as 
shown in Fig. [TJ 




Figure 2: (Color online) The transition pressure Pt as a functionof L and pt by using the 
dynamical approach with and without the parabolic approximation in the MDI interaction. 
Taken from Ref. [32]. 



The transition pressure P t at the inner edge of the neutron star crust is also an 
important quantity that might be measurable indirectly from observations of pulsar 
glitches [U [TTJ. Shown in Fig. [2] is the P t as a function of L and p t by using the 
dynamical approach with both the full MDI EOS and its PA. Again, it is seen that 
the PA leads to huge errors for large (small) L (p t ) values. For the full MDI EOS, the 
P t decreases (increases) with increasing L (p t ) while it displays a complex relation 
with L or p t in PA. From the constrained L values, the value of P t is limited between 
0.01 MeV/fm 3 and 0.26 MeV/fm 3 . 
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4 Improved constraint on the radius-mass relation 
of neutron stars 



The constrained values of pt and Pt have important implications in many properties 
of neutron stars [H [8] [12], [27]. As an example, we have examined their effect on 
constraining the mass-radius (M-R) correlation of neutron stars. The crustal fraction 
of the total moment of inertia AI/I of a neutron star can be well approximated 

AI _ 28irP t R 3 (1 - 1.67f - 0.6£ 2 ) 
T~ ~ 3Mc 2 \ 

where m\> is the baryon mass and £ = GM/Rc 2 with G being the gravitational 
constant. As stressed in Ref. [3], AI/I depends sensitively on the symmetry energy 
at subsaturation densities through P t and pt, but there is no explicit dependence on 
the higher-density EOS. So far, the only known limit of AI/I > 0.014 was extracted 
from studying the glitches of Vela pulsar [11] . This together with the upper bounds on 
Pt and pt (pt = 0.065 fm~ 3 and P t = 0.26 MeV/fm 3 ) sets approximately a minimum 
radius of R > 4.7 + 4.0M/M Q km for the Vela pulsar. The radius of Vela pulsar 
is predicted to exceed 10.5 km should it have a mass of 1.4M Q . We notice that a 
constraint of R > 3.6 + 3.9M/M Q km for this pulsar has previously been derived in 
Ref. [IT] by using p t = 0.075 fm -3 and P t = 0.65 MeV/fm 3 . However, the constraint 
obtained in our study using for the first time data from both the terrestrial laboratory 
experiments and astrophysical observations is more stringent. 

The above constraints are shown in Fig. [3] together with the M— R relation ob- 
tained by solving the Tolman-Oppenheimer-Volkoff (TOV) equation. In the latter, 
we have used the well-known BPS EOS [6] for the outer crust. In the inner crust with 
Pout < P < Pt, the EOS is largely uncertain and following Ref. [26], we use an EOS of 
the form P = a + foe 4 / 3 with the constants a and b determined by the total pressure 
P and total energy density e at p out and p t . The full MDI EOS and its parabolic 
approximation with x = and x — — 1 are used for the uniform liquid core with 
p > Pt- Assuming that the core consists of only the npe matter without possible new 
degrees of freedom or phase transitions at high densities, the PA leads to a larger 
radius for a fixed mass compared to the full MDI EOS. Furthermore, using the full 
MDI EOS with x = and x = — 1 constrained by the heavy- ion reaction experiments, 
the radius of a canonical neutron star of 1.4M Q is tightly constrained within 11.9 km 
to 13.2 km. 



1 + 



2P t (l + 5£-14e 



p t m h c 2 i 2 



(5) 
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Figure 3: (Color online) The mass-radius relation M-R of static neutron stars from 
the full EOS and its parabolic approximation in the MDI interaction with x = 
and x = —1. For the Vela pulsar, the constraint of AI/I > 0.014 limits the allowed 
masses and radii. See text for details. Taken from Ref. [32] with small modifications. 



5 Summery 

Using the MDI EOS of neutron-rich nuclear matter constrained by recent isospin 
diffusion data from heavy-ion reactions in the same sub-saturation density range as 
the neutron star crust, we have determined the density and pressure at the inner edge, 
that separates the liquid core from the solid crust of neutron stars, to be 0.040 fm -3 
< Pt < 0.065 fm~ 3 and 0.01 MeV/fm 3 < P t < 0.26 MeV/fm 3 , respectively. These 
constraints have allowed us to determine an improved mass-radius relation for neutron 
stars. Furthermore, we have found that the widely used parabolic approximation 
to the EOS of asymmetric nuclear matter leads to significantly higher core-crust 
transition densities and pressures, especially for stiff nuclear symmetry energies. 
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